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Abstract 

We construct a 3+1 dimensional holographic model dual to a parity violating 
hydrodynamic system in 2+1 dimensions. Our model contains gravitational and 
electro dynamic Chern-Simons terms coupled to a neutral pseudo scalar 6, and a po- 
tential composed of quadratic and quartic terms in 6. The background is a charged 
black brane. We study the hydrodynamics to first order in spacetime derivatives 
near the probe limit of the pseudo scalar, by extracting the transport coefficients 
from the scalar, vector, and tensor modes of bulk perturbations. We study two 
mechanisms for breaking the parity of the boundary fluid: the parity is either spon- 
taneously broken by the nonzero vev of the dual pseudo scalar operator, or by the 
pseudo scalar source on the boundary. The resulting transport coefficients have 
unusual behaviors compared to those in ordinary materials. 
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1 Introduction 



In the last decade, AdS/CFT correspondence has been demonstrated to be a useful tool 
to study the strongly coupled gauge theories. One of the particular applications recently 
attracting significant attentions is holographic hydrodynamics [UI2]- The hydrodynamics 
that is particularly well described by this approach so far is relativistic in nature. By 
applying this technique, various hydrodynamics transport coefficients of certain classes of 
strongly coupled conformal field theories have been computed subject to the underlying 
symmetry. According to this holography, the near horizon physics of the regular black 
hole in an asymptotically AdS spacetime is able to describe the long wavelength dynamics 
of a strongly coupled finite temperature field theory on the boundary. 

It is well known that the conventional hydrodynamics can be constructed from the 
underlying symmetry of the system without knowing much details about the microscopic 
degrees of freedom. In the hydrodynamic regime, a system can be completely described by 
a set of conservation equations, the relativistic generalization of the well-known Navier- 
Stokes equations. If the system has some dissipative parts in the conservation equations, 
their form are subject to the second law of thermodynamics. One interesting feature 
of the aforementioned holographic method [2] is that it turns out to incorporate the 
constraint of positivity of entropy production, and successfully reproduces the relativistic 
hydrodynamic equations without referring to this constraint. It has been conjectured by 
an explicit example that the solutions of the relativistic hydrodynamic equations are in 
fact in one to one correspondence with the regular long wavelength solutions of the bulk 
Einstein equations in AdS spacetime. Afterwards, the holographic technique had been 
applied to a wide range of strongly coupled systems dual to various kinds of asymptotic 
AdS black hole backgrounds [31 IU El [6] , and further generalized to describe the superfluid 
phase [7|. 

Besides being an elegant complimentary approach towards deriving the hydrodynamic 
equations, the holographic approach provides us new transport coefficients which do not 
appear in the conventional phenomenological construction. This essentially inspires us 
to reconsider the whole approach of constructing the hydrodynamics from more general 
perspective. Positivity of the local entropy production is an important guiding principle 
for specifying the constitutive relations among the thermodynamic variables. It is first 
pointed out in (8] that, in order to accommodate those new terms in the stress-energy 
tensor and the global current, the entropy current has to be modified. Based on this, 
attempts have been made to construct the most general relativistic hydrodynamics order 
by order in derivative expansions, to incorporate various field theory anomaly- induced 
transport coefficients O [10] . See also [TT] for the recent proposals of constructing hydro- 
dynamics from the field theoretic perspective without referring to the entropy current. 

The holographic approach can be applied to the systems with parity violation [T21 EH 
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Q31 [T5]. Recently, parity violating transport in high energy theory has gained significant 
interests, due to the charge asymmetry observed in the heavy-ion collision experiments 
such as STAR [E], PHENIX pi] and ALICE |18j. The parity violating chiral magnetic 
effect (CME) was proposed to explain the aforementioned observation and subsequently 
under intensive study during the past few years [T§|. On the other hand, the low energy 
transport properties of parity violating systems had long been the subject of interest in the 
condensed matter research. Anomalous Hall conductivity is one of such effects which had 
been measured experimentally, see [20] for recent review. The thermal Hall conductivity 
in the presence of background magnetic fields is demonstrated to be a useful probe for 
the quasiparticle transport in the high temperature superconductors [21] • There are also 
discussions on the anomaly- induced thermal Hall effect in topological insulators [22]. A 
unified field theoretic model had been proposed for the anomalous electric, thermoelectric, 
and thermal Hall effects in the ferromagnetic materials [23]. In hydrodynamics, the Hall 
viscosity, which is the parity violating analog of shear viscosity, is already studied in the 
field theoretic approach [2H 125] . 

When the parity of the underlying theory of a system is broken, extra transport 
coefficients arise in the long wavelength limit. For the strongly interacting fluids, these 
new coefficients are difficult to evaluate from field theory perspective. One well-studied 
example is the Hall conductivity. From the field theory point of view, the Hall conductivity 
may originate from a bare Chern-Simons term in the fundamental theory, or the Chern- 
Simons term induced by radiative corrections. Even in the lattice computation which is 
very useful when the coupling is strong, there is a sign problem occurring to the Chern- 
Simons term after performing the Wick rotation. However, the AdS / CFT correspondence 
offers a powerful non-perturbative tool to calculate these new transport coefficients. In 
this paper, we use this method to compute all possible transport coefficients in first order 
derivative expansion of a holographic, strongly coupled, charged fluid in 2 + 1 dimensions 
with parity violation in the hydrodynamic limit. 

This paper generalizes our previous construction in [13] by considering Chern-Simons 
terms coming from the electromagnetic sector, FF, and the gravitational sector, RR, both 
coupled to a neutral pseudo scalar field. Our four dimensional gravitational background is 
a charged black brane in the probe limit of the pseudo scalar. The parity of the boundary 
fluid is broken by the operator dual to the bulk pseudo scalar field. Unlike [H] which 
studies the holographic parity-violating hydrodynamics by means of the linear response 
functions, we follow the standard fluid-gravity correspondence technique [2] to compute 
the hydrodynamic transport coefficients, including those parity-odd ones, and found some 
exotic properties of the latter which are yet to be understood from field theory perspective. 

We employ two different parity breaking mechanisms in our system, corresponding 
to two different boundary conditions for the bulk pseudo scalar field. In one scenario, 
the pseudo scalar is sourceless on the boundary, and the parity is broken spontaneously 



2 



by the nontrivial vacuum expectation value of the dual scalar operator developed below 
certain critical temperature T c . As a result, all parity violating coefficients are related to 
the order parameters for the second order phase transition of our system, and becomes 
non-zero below T c . In this case, the boundary theory is still conformally invariant, leaving 
the bulk viscosity zero. Our result shows that, some of the parity-odd coefficients, such 
as the anomalous electrical Hall conductivity a and Hall viscosity ?}a, go to zero with the 
mean field critical exponent (T c — T) 1 / 2 , while others, including the anomalous thermal 
Hall conductivity k and the anomalous heat Hall conductivity £, diverge with the critical 
exponent (T c — T)~ 3 / 2 at T c . These properties are novel, and it would be interesting if 
they can be realized by field theory models. 

In another scenario, the pseudo scalar is sourced on the boundary, and the parity of 
the hydrodynamic system is broken by the externally applied source of the dual boundary 
pseudo scalar operator. Moreover, the source also breaks the conformal symmetry. There- 
fore, in the hydrodynamic limit, in addition to the parity violating coefficients in the first 
scenario, new coefficients such as the bulk viscosity £, curl viscosity (a and magnetic vis- 
cosity £b also arise. We explore the holographic realization of parity breaking with both 
scenarios in detail, and study the finite temperature behaviors of various hydrodynamic 
transport coefficients of a strongly coupled system. 

The structure of this paper is organized as follows: in section [2] we briefly review 
the 2 + 1 dimensional parity violating hydrodynamics, mostly following Ref. [T3j. The 
notations subsequently used in our paper are also introduced here. In section [3] we present 
our holographic model and the ansatz for the derivative expansions of perturbations. The 
background solution, including two different boundary conditions for the pseudo scalar 
corresponding to two distinct mechanisms of breaking the boundary parity, is given in 
Section HI We further reproduce and analyze the transport coefficients of parity violating 
hydrodynamics from the bulk gravitational theory, under the probe limit of the pseudo 
scalar field. Our research is concluded in Section The first order perturbation equations 
of motion of the bulk theory before taking the probe limit are listed in the Appendices. 



2 Hydrodynamics with broken parity in 2 + 1 dimen- 
sions 

In the hydrodynamic limit which assumes local thermodynamical equilibrium, a system 
with global U(l) symmetry under a non-dynamical external gauge field strength F^f 
satisfies the following local conservation equations: 

V M T"" = FZJ» V = 0, (1) 

without referring to the microscopic degrees of freedom, where T^ v and J M are the energy- 
momentum tensor and the conserved U(l) current respectively. In a 2+1 dimensional 
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parity violating hydro dynamic system described by the velocity u p , temperature T, and 
chemical potential /i, the most general constitute equations of T pu and J p in principle can 
be written down in terms of all possible combinations of these local macroscopic variables 
which respect the underlying symmetry and the constraint of positive entropy production 



= Corfu" + P A^ - (CV a m a + Ca^ + (bB)A^ - r)a^ - ?j A a^, (2) 

,^ + ae^u u E p + K6^u v V p ^ 



J» = pu » + a E»-K^V u ^ + ae» up u y E p + Ke^ p u u V p ^ + ^ vp u u V p T. (3) 



where the combinations of the local variables are 



Q = -e^Pu^VuUp (vorticity) 



i 1 

B = —~^ up u p F^ x p t (external magnetic field) 



E p = u\F P xt (external electric field) (4) 

A pu = if + u p u u (projector to the direction transverse to u p ) 

= A" a A^(V^ + V^-f, a/j V/) 



Note that there is an ambiguity in interpreting the transport coefficients in the constitution 
equations unless a particular choice of out-of-equilibrium definition of energy density, 
charge density, and local velocity is given, i.e. a frame is chosen. In our paper we choose 
the Landau fram^_, such that u l is identified to be the energy transport velocity T OT , i.e. 
T pv u v = —eou p and J^u^ = —p. Then eo, p, Po are interpreted as the energy density, 
charge density, and pressure respectively in the equilibrium configuration where B = 0, 
Q = 0. The local velocity is normalized as u p u p = — 1. 

The coefficients in the constitutive equations are interpreted as follows. £ and r) are 
the bulk and shear viscosities, a is the electric conductivity; k = oT is the thermal 
conductivity. All these above are the dissipative transport coefficients appeared in the 
usual parity-even system, and they must be positive. The remaining coefficients (with 
tilde) arise from the parity violating effect. They are dissipationless and there is no 
constraint on the signs of their values. Among them, (a and ?}a are the curl and Hall 
viscosities respectively. A cartoon depicting the physical pictures of (a and tja is given 
in Fig. 1 of 113]. In this paper we call (b the magnetic viscosity, a is the anomalous Hall 
conductivity^ which describes the Hall effect without the presence of external magnetic 
field. We name k and £ the anomalous thermal Hall conductivity and anomalous heat 



1 Here the indices {/i, v\ run through all indices of the 2+1 dimensional spacetime, while {i,j} only 
label the spatial dimensions. 

2 We use different but equivalent parametrization in ([3]) compared to that in [TJ]. Our a is in fact 
a + xe in [14] . Their \e is equal to our a — y. 
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Hall conductivity respectivel}|§, which describe the transverse current due to gradients of 
^ and T without external magnetic field. 



[T4"j refers to £, 77, a, tja and a as transport coefficients, and to (a, (b, k and £ as ther- 
modynamic response parameters. The difference is, in the Kubo's formula formulation, 
those thermodynamic response parameters are computed by setting zero-frequency before 
taking zero-momentum limit, while for those transport coefficients the order is reversed. 
With these details in mind, for convenience we just call all of them transport coefficients 
in the rest of our paper as the convention does. 



3 The holographic set-up and the derivative expan- 
sion 

Our bulk model dual to the 2+1 dimensional parity-odd hydrodynamics is described by a 
3+1 dimensional Einstein-Maxwell-pseudo scalar Lagrangian with a negative cosmo logical 
constant. The U(l) gauge field is sourced by the gravitational background, and the neutral 
pseudo scalar couples to the gravitational and electrodynamical Chern-Simons termj^: 



where 



PR _ pM PQ pN pM PQ ._ 1 C PQRS pM 

titi — K n ti MPQ , ti n ■— — £ ti NRS-, 



FF = F MN F MN , F MN :=~e MNP QF PQ , 

with e MNp Q being the usual 4-dimensional Levi-Civita tensor. We assume that the gravi- 
tational and the electrodynamical Chern-Simons coupling constants are of the same order, 
in order for the parity breaking effects arising from both terms to appear at the same order 
of the probe limit we take later. For simplicity, we take A = A' in (jSJ). 

An appropriate boundary action S^dy can be added to the bulk action to make the 
variational principle consistent on the boundary: 

S Uy = 2 / d 3 xVh ( — - \ 9 n A e ABCD K E B W c K DE \ . (6) 
J \lbirG N 2 J 

We define as the outward pointing normal vector to the asymptotic boundary of 
AdS with unit magnitude. Kab = hA C hB D V c^d is the extrinsic curvature tensor of 
the boundary slice, with the induced metric h A B = 9ab — katie, and K = H ab Kab is 



3 In this paper, the term "anomalous" means that the response is transverse to the input without 
the presence of background magnetic field, e.g. anomalous Hall effect. On the other hand, we use 
"anomaly-induced" to denote the parity-odd effect due to the gauge anomaly, as in Section [1] and [5] 

4 the bulk electromagnetic coupling e 2 is set to one. 
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the trace of Kab- According to the prescription of the AdS/CFT correspondence, the 
on-shell variation of the boundary action with respect to the induced metric gives rise 
to the energy-momentum tensor of the dual field theory. Additional boundary terms 
involving 9 originated from 9RR in the bulk action in principle modifies the boundary 
energy-momentum tensor [T21I2B] . But as long as 9 is a relevant perturbation, i.e. m 2 < 
such that 9 is asymptotically normalizable, these additional boundary terms vanish and 
have no contribution to the hydrodynamic transport coefficient arising from the energy- 
momentum tensor. See the Appendix of [12] for detail. In this paper we focus on m 2 < 0, 
and set 1/16ttGn = 1 in the bulk action. 

The 3+1 dimensional bulk coordinates are denoted by x M = (r, x^), where x^ = 
(v, x, y) are the coordinates labeling the boundary slice. 

The parity in the bulk Lagrangian ([5]) is conserved. Parity violation in the boundary 
theory can arise either spontaneously from the 9 condensate corresponding to the bulk 
solutions which satisfy the sourceless boundary condition, or "explicitly" from the bulk 
pseudo scalar with a source on the boundary. The potential V(9) in our model takes the 
fornffl 

V(p) = ±m 2 P + \cP. (7) 
The equations of motion arising from ([5]) and ([7]) are 

Rmn — ~9mnR + ^9mn — ^Cmn = Tmn{Q) + Tmn(A), 

V M F MN = Xd M 9F MN , (8) 



where 

T MN (0) = \d M 9d N 9-^g MN {d9f -^g MN V{9), 

Tmn(A) = -Fm A F~na — -9mnFabF ab '', 

C MN = (Wp9)e PQ ^ M \/sR N) Q + (W P \/ Q 9)R^ MN ^ p . 

Cmn is a symmetric traceless tensor analogous to the Cotton tensor in three dimensions, 
and the symmetrization is defined as RQ( MN ) p ■= \{EP MNP + EP NMP ). 



5 It is known that with the potential (T7J), a neutral hairy scalar can only satisfy the sourceless boundary 
condition in a chargeless AdS-Schwarzschild background for certain range of c in which the positive 
energy theorem is violated [2D] ■ Introduction of the gauge fields sourced by the gravitational background 
eliminates this problem. 

6 In this paper, for the case where 9 is asymptotically sourceless, we focus on the regime — | < m? < — | , 
as in order for 9 to have a non- vanishing vev on the boundary, it has to violate the Breitenlohner-Frecdman 
bound in the near horizon region, which is m 2 = — | for our black brane background whose near horizon 
geometry at T = is AdS2 ■ For the case where 9 is sourced on the boundary, we take the m 2 range such 
that 9 is dual to the relevant deformation in the boundary theory, i.e. — | < m 2 < 0. 
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The background solution to the bulk e.o.m.()H]), which describes a uniform ideal fluid 
on the boundary with velocity w M , temperature T and chemical potential /x, is given by 
the following boosted charged brane ansatz 

ds 2 = -2H(r,M,Q) Ufl dx^dr - r 2 f(r, M,Q)u fl u l/ dx (1 dx u + r 2 A^dx^dx u := ds 2{0 \ 
9 = 6(r,M,Q) := 9®, (9) 
A = [A{r,M,Q) Ufl + Al xt ]dx^ := A®, 

A ext ( a ext a ext \ext\ 

A IJ, — K^-v >Ac i^-y J' 

where m m is the three- velocity on the boundary with = — 1, fi is given by the difference 
between the boundary and the horizon value of A v , T is the black hole temperature, and 
A M „ = Tj^ v + u^Uy is the projector to the directions transverse to w M . Here we also turn 
on the constant external gauge fields A ext , besides the dynamical ones. Since A ext are 
constants in our setup, they have no contribution to the background bulk dynamics. 
However, the external electric and magnetic fields will arise from perturbing A ext , which 
will be clear shortly. The superscript ^ denotes the background ansatz with respect 
to the perturbations and corrections introduced later. (See ([T5]) .) We also choose the 
Landau frame for our bulk background. 

For 2+1 dimensional point of view, by promoting the uniform local variables w M , T and 
fi to slow-varying functions of spacetime without breaking local conservation laws, one can 
add dissipative and non-dissipative correction terms with various transport coefficients 
to the energy-momentum tensor and the current of the uniform ideal hydrodynamics. 
One can also turn on external electric and magnetic fields which also contribute to the 
corrections to T^ u and J^. 

In terms of the dual gravity, this is equivalent to promoting the constant bulk quanti- 
ties M, Q, uF" and A ext to slow- varying functions of the boundary coordinates x M , and then 
determining their effective dynamics. Due to this field promotion, the original bulk fields 
undergo derivative expansions, which together with the correction ansatz solve the e.o.m's 
order by order. By solving all the equations of motion at each order, the boundary hy- 
drodynamics in terms of these slow- varying fields are realized by the standard AdS / CFT 
correspondence, and the transport coefficients are extracted from the asymptotic behav- 
iors of the bulk field solutions[2]. For the charged fluid case, the holographic description 
is summarized in [B]. We will apply this approach to our bulk theory (jSj) to study the 
dual parity violating hydrodynamics on the boundary. 

In the co-moving frame, at the origin of the boundary coordinates x M = we have = 
(1,0,0). The promoted slow varying fields u' 1 (x u ), M(x u ), Q(x u ), A ext (x u ) are expressed 
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in terms of the Taylor expansions with respect to the origin: 



u»{x v ) = (1, x»d u /3 l ) = (1, 5/3*), 

M(x u ) = M + x v d v M, (10) 

Q{x v ) = Q + x u d u Q, 

A e *\x u ) = Af + x v d v Af. 

As a result, the bulk metric components f(r,M,Q),H(r,M,Q), dynamical gauge fields 
A(r, M,Q), and pseudo scalar 9(r,M,Q) also receive derivative expansions, 

6F{r, M,Q) = ^ x^M + |^ x»d,Q, (11) 

where F = {/, if, A, 8}. Since such perturbed ansatz is no longer a solution to the original 
equations of motion, the following corrections to the bulk fields must be introduced to 
compensate the deviation from the original solution due to derivative expansions f llUp and 

ds 2 (1) = r 2 k(r) dv 2 + 2Hh(r) dvdr + 2r 2 j i (r) dvdx 1 - r 2 h(r) dx % dx % + r 2 otij(r) dx l dx j , 
A« = a v {r) dv + ai{r) dx\ (12) 
0« = p(r), 

where CKy is symmetric traceless. 

Combining (1§|)~( 1T2"|) . the overall charged black brane ansatz up to first order pertur- 
bation reads 

ds 2 = ds 2 ^ 

+ e [-r 2 5f dv 2 + 2 5H dvdr - 2 r 2 (l - f(r)) 5/3* dvdx 1 - 2 H(r) 6/3* drdx*] 
+ eds 2{1 \ (13) 

A = A (0) + e (-6Adv + A(r)5(3*dx* + A (1) ), 

6 = 6 {0) +e(56 + 6 {l) ). 

The parameter e is introduced to keep track of the orders of derivatives and corrections. 
The superscripts and W in and (fT2"j) label the 0(e°) and 0(e x ) ansatz respectively, 
and the whole ansatz (fl3|) solve the bulk e.o.m.'s up to 0(e l ). One can study the higher 
order hydrodynamics by performing such procedures order by order, to solve the 0(e n+l ) 
equations of motion with 0(e n+1 ) ansatz provided the 0(e n ) equations are solved by the 
0(e n ) ansatz. 

In the next section, we will present the solutions at the probe limit of the pseudo 
scalar field, and extract the corresponding transport coefficients for the boundary fluid. 



S 



4 Transport coefficients of parity violating fluids in 
2 + 1 dimensions from gravity 

In this section, we will solve the bulk equations of motion up to 0(e) perturbations, and 
extract the first order transport coefficients of the dual hydrodynamics. 

It is clear from (JS} that the gravitational background receives backreaction from non- 
trivial pseudo scalar profile. Despite that the full analytic solution with backreaction is 
not yet known, we can still explore the dual hydrodynamics in the probe limit of the 
pseudo scalar field 

A0 , V(9)^XV(9), A^O. (14) 

We will see that, in vector and tensor modes, by taking the probe limit as in (]14p . we are 
able to separate parity-even transport coefficients which can be see at 0(A°), from the 
parity-odd transport ones which only appear at higher order in A. 

We first present the background solution in the probe limit. In the Einstein and 
Maxwell equations, since the leading order term are of 0(A°) while the 9 terms contribute 
at 0(A 2 ), the pseudo scalar decouples at the leading order of the A expansion. The 
background is therefore a charged black brane in AdS^, 

ds 2m = 2dvdr-r 2 f(r)dv 2 + r 2 (dx 2 + dy 2 ), 
M O 2 

A (o,6) = JQ d 

r 

where another superscript ^ is employed to label the order of A, while the untilded one 
( n ) denote the order of e. The horizon is at r — ru, ancj^l M, Q are the black brane mass 
and electric charge respectively. The metric is asymptotically AdS^ with the AdS radius 
taken to be 1. The temperature and the chemical potential are 



3r H / Q 2 \ 2Q 

T = 17( 1 -3^J- » = - H (16) 



The probe limit (fl"4"l) implies that the scalar field is in fact 9 = #( 0,1 ). Since ffl5|) gives 



7 The action and equations of motion have the following scaling symmetry: 

r br, (v,x,y) - (v,x,y), 
Q -4 b 2 Q, M b 3 M, 

One can make use of this invariance to rescale the horizon to Th = 1 by setting b = l/rn, M = 1 + 'in 
and Q 2 = 3k, but the independence between temperature T and the chemical potential fi (or the charge 
Q) in the boundary theory will be lost. 
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rise to vanishing RR and FF terms, the Klein-Gordon equation in (jSJ) reduces to 

r + { f + t)f-^ = . (17) 
/ r r 2 f 

Since f(r = th) = 0, it follows straightforwardly from the above equation that on the 
horizon 9 is subject to the condition 

dV/dO 



9'( 



r 2f, 



'If 



On the other hand, the pseudo scalar behaves asymptotically as 



J(°> (O) (0) . 3 J [9 



= + + •••, A± = - ± W- +m 2 . (19) 

For — | < < 0, is identified as the boundary source of 9, while (C) (0) is the response 
(or, condensate) on the boundary dual to the source J(°\ For —9/4 < m 2 < —5/4, one is 
free to choose either or (O) 1 - ^ as the source and the other as condensate [21] • In this 
paper we choose as the source and (O)^ as the condensate. 

In our model, we consider two types of asymptotic boundary conditions for 9, corre- 
sponding to two mechanisms to break the parity in the boundary theory: 

(i) Boundary condition with a source: 

One can turn on the boundary source for 9, such that the condensate {O}^ is 
determined by the regularity condition on the horizon. In the analysis throughout this 
paper, We set = 1 for convenience. In this case, the boundary parity and conformal 
symmetry are broken by the source 

(ii) Sourceless boundary condition: 
The other scenario is, the source 

J(0) 

can be switched off consistently, such that the 
parity and conformal symmetry in the boundary theory are broken spontaneously by the 
nontrivial vev of (O)^ developed below certain critical temperature depending on m 2 . In 
this case, from the Klein-Gordon equation (|17[) . one can derive the mean field description 
of (C) (0) ~ (Tc-T) - 5 near T c . This is confirmed by the numerical results^! in Fig. [TJ 

Having the unperturbed background solution in hand, in the following we solve the 
0(e) equations of motion, in terms of the scalar, vector, tensor modes with respect to the 
boundary spatial SO (2) symmetry in the probe limit. The full 0(e) equations with back- 
reaction are presented in the Appendices for interested readers. We will denote the 0(e) 
Einstein, Maxwell and Klein-Gordon equations by E^j N , and KG^ respectively. 



4.1 Scalar mode 

In the probe limit (Till) and tp — > \<p, there are four dynamical equations for the scalar 
modes arising from E$, H(r)M v ( 1 \ ^J-E„ + E$ and KG^\ The first three equations 

8 We take c = 0.5 in the potential V{9) throughout the numerical analysis in this paper. 
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Figure 1: - as a function of T//x with linear scale (on the left) and log scale (on 

the right) for various m 2 . In these figures 9 is asymptotically sourceless, i.e. = in 

(]19[) . - vanishes at and above the critical temperature, with a critical exponent 

equal to 0.5, as being indicated by the right plot. 



are 

1 



1 

r 



2K r 2 a v ')' + 2A'ti = 0(A 2 ), (20) 



^ k y + M {r 'fl' h )'-l A ' a °' = ;^ + o(A a ). 

Substituting f(r) and A(r) from O(e ) solution ( I15p into the equations above, one obtains 
h(r) = Cl -^ + 0(A 2 ), 

<r) = c 3 -^-^ + 0(A 2 ), (21) 

U \ O , °5 c 4 Q c 2 M c 2 Q 2 1 „ , 2 . 

fc(r) = -2ci + — + — - — — + — — + -diPt + 0(A ), 

where ci, . . . , C5 are arbitrary constants. 

In our model, the integration constants C\,cz are set to zero by choosing the normal- 
izable boundary conditions for h(r) and a v (r). 04,05 takes zero value due to the fact we 
choose the Landau frame such that tt^ux = and v^u^ = 0, where tt^, z/ M are O(e) 
corrections to the energy momentum tensor T mA and the current respectively [14] . c 2 
can be absorbed by redefinition of the coordinate r and set to zero [2J |5]. Hence, 

h{r) = 0, a v (r) = 0, k{r) = - (22) 

r 

With (J22D, the constraint equations (jTSJ) and QBE) arising from g rr E$ + g vr E$ and M r ^ 
reduce to 

d v Q = -Qdifa + 0(A 2 ), (23) 
«9„M = J-M-dA + 0(\ 2 ). 
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By substituting in (j22p and (I23p . the fourth dynamical equation, i.e. KG^ of 0(A), reads 



d 2 V 



dd' 



-tp 



1 , 2n,M 3M / dQ V 2Q ( dQ 
- (r z Q) 1 



r V dM J r \ dQ 



+ 



2,4,4' 



24' 



5, (24) 



where 5 = = d x A e y xt - d y A e x xt . Note that unlike ([21]), the KG® equation (1231) 

describes parity violating effects. 

To solve i^G^ 1 ), we decompose y> into the "bulk", "curl", and 'magnetic" parts which 
are proportional to d^, di(3j, B respectively, as they are independent from each other: 



v( r ) = M r ) diPi + ¥c{r) (d x (3 y - d y (3 x ) + ip B (r) B. 



(25) 



such that 



dO 2 
d 2 V 



1 2 , 3M / dQ V 20/ <9# \ 



r \ <9M 
2AA' 



^ 4 M)'- d6)2 

p^M*) dff 2 = 
(fib, tp c and ifB are subject to the regularity condition on the horizon 

, 1 



r dQ) 



(26) 



2A' 



{r H ) 



[r H ) 



r 2 fl 
1 

1 

^77 



1 d , 3M d / dQ \ 2Q d ( dQ\ d 2 V 

r dr 

1 d 

4r 2 dr 



r + 

r dr \ dM J r dr \ dQ J dd 2 



, _~ 2AA' d 2 y 



<Pb 



r=rn 



dQ 2 



(27) 



r=r H 



2 A' d 2 V 



+ 



dQ 2 



r=r H 



and asymptotically it behaves as 



J (O) 



J(°) + e( J b (1) ^ft + J c ( % d^ + J^B) 



(O) (0) + e{(0)? } m + <W Wj + (C)^ 5) 



(i) 



1 (O) (0) + e«0) ^ + <0> ^ m + (<D)$ B 



A 1 ) 



(1) 



(2« 



where = and we have chosen the boundary condition { = 1, = 0} such 
that {(O)^ , (O) 1 ' 1 ^} are determined by the regularity condition of Q and ip on the horizon 
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respectively. The boundary conformal symmetry and the parity are broken by the source 
j(o) c hoose = 0, because the source on the the boundary is specified by 0(e°) 
background and we don't want to be modified by the bulk perturbation^. Then, the 
bulk viscosity (, curl viscosity Ca, and magnetic viscosity Cb can be extracted from the 
following expression of the trace of the boundary energy-momentum tensor 



(T^) (1) = A„(J(°)(0) (1) + JW(0) (O) ) 

= A^^O^dA+iO^e^d^ + iO^B) (29) 
= 2Cd i /3 i + 2CAe ij d i /3 i + 2C B B, 



where 



c = 


~Y 


«> 


Ca = 


A_ 

~Y 




Cb = 


A_ 
~2~ 





(bulk viscosity) 

(curl viscosity) (30) 
(magnetic viscosity) 

The numerical results for the dimensionless combination -, — and are in Fig. 

S ' S p/^l ° 

EHU Our numerical analysis shows that (/s is positive, diverges at T//i = and tends 
to a certain finite number at large T/fi. Ca/ s is negative and asymptotes to a certain 
finite value at large T/fi. It becomes more negative toward smaller T j \x before it starts 
to increase and eventually flips sign at very small T//i. The T/fi point at which Ca/ s 
vanishes in Fig. [3] turns out to be the same for different m 2 . This property is novel and 
the underlying physics is unclear. However, since the curl viscosity is dissipativeless, the 
change of sign is allowed, as there is no constraint on the sign of Ca- Ca/ s approaches 
certain finite value when T/fi = 0. On the other hand, ^ is negative and decreases 
monotonically as T/fi reduces, until it hits certain finite value at T//i = 0. 



4.2 Vector mode 
4.2.1 Parity-even part 

For the O(e) vector mode perturbations, there are two dynamical equations from E ri ^ 
and r 2 M %< ^\ In the probe limit they reduce to 

^faf^-r'A'ji = -^Q-A'd^ + OiX 2 ), (32) 

9 On the other hand, if one chooses both JW = and — 0, the viscosities arising from the scalar 
mode are just trivially zero in (f29)) . 



13 




Figure 2: (/s versus T/fi for various m 2 . The right figure is for the very small T/fi. The 
boundary conditions of 9 and ip underlying these numerical results are summarized below 
equation (128]) . 




Figure 3: (a/s versus T//x for various m 2 . The right figure is for the T//x range close 
to zero. The boundary conditions of 9 and ip underlying these numerical results are 
summarized below equation ( I2"g|) . 




Figure 4: jj^ versus T//i for various m 2 . The right figure is for very small T/jj. The 
boundary conditions of 9 and ip underlying these numerical results are summarized below 
equation (|28|) . 



14 



where r 2 ji = g$. The dynamical equations (l3Tj) and (132]) describe a parity-even system. 
In the probe limit, the constraint equation from ^Fyp + E$ reads 

diM = -3Md v f3 i + 2QE l + 0(\ 2 ), (33) 

where E i = F^f 1 is the external electric field. 

a, and jj can be solved analytically from (|3~Tj) and ( 13"2"j) . After integrating both equa- 
tions once, one can write down the decoupled differential equation for jV 



(r 6 /)if + (4r 5 /)j- - 4Q 2 J, = (-2r 3 / + ^)<%A + + 2Qb, 



(34) 



where b is an undetermined integral constant. Observing that the left hand side can be 
written into a total derivative r 2 (r 4 f 2 (ji/f)')', one obtains the solution for ji, 



Ji 



f / ds 



1 



s'f 2 



(Q 2 2M 2 \ '*4 * n 



2Q 



2Qb\ 
s 5 f 2 ) 



(35) 



where we have integrated from oo to r to derive the solution ji from f l34j) [35J. At the 
asymptotic infinity where / — > 1, the closed form for jj expression can be written down 
(where the % index on the right hand side is supressed): 



9«A 



+ o 



(r 5 )' 



(36) 



where 64 = ^- d v (3i + -Q. This is consistent with the normalizable boundary condition 
which demands vanishing r° coefficient, and the requirement of the Landau frame choice 
that r~ 3 coefficient (i.e. the boundary expectation value corresponding to ji) has to 
vanish. Finally, b is fixed by the regularity condition on the horizon, 



and ji{r — th) is 



3i{rH) 



1 



2Q 2 1 



(4Qd v Pi + 2diQ), 



1 



3M r H 



3M r 2 H 



r H v 3M r H J 

Next, di(r — > 00) is obtained from integrating equation (13T1) from infinity to r: 



(37) 



(3J 



«i(r^ 00) = ^£^{(^)' + 2^A} = ^[3 4 j;' 



r 
«Joo 



1 

2Q 



1 2 a 



(39) 



With the standard AdS/CFT correspondence, one can derive the boundary current 
from the bulk solution by 



hm — - 

r->oo 



A" 

4Q 2 1 
3M 



lim 

r— >oo 



F* + 



2 



E l + -b, 



2Q 2 1 



2Q 

3Mr H ) \3M 



diM - diQ 



(40) 
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2 4 6 8 10 12 fi 0.1 0.2 0.3 0.4 0.5 0.6 



Figure 5: a (the top figures) and -^h (the bottom ones) versus ^, for the ranges of high 
— (on the left) and low — (on the right), k and a are related by k = crT. The value of a 
and -j- depend only on the charged black brane background but not on 9. 



This equation exactly corresponds to the current of the fluid 



J^alE.-Td^), (41) 



T 

if we identify the temperature T, chemical potential /i as in (1161) . and the electric conduc 
tivity a, thermal conductivity n as 

4Q 2 1 V _ fAnr 2 H T 



3Mr H J V 3M 
1 / 4Q 2 1 \ 2 f3M 4Q'< 



K - aT --A l -^- H ) (42) 

The thermodynamic relation of the current for the fluid in (|41|) is reproduced from gravity. 
Note that in a, T is also a function of r# and M, and it is straight forward to conclude 
from the above expression that a — > 1 as T — > oo and a — > as T — > 0. The temperature 
and r# dependence in a is exactly the same as that derived in [34] by Kubo's formula. 

The numerical behaviors of dimensionless a and -A-^ as functions of T/fi are displayed 
in Fig. [5j This figure shows that under fixed chemical potential, the electric conductivity 
decreases as the temperature is reduced which is different from the usual metal property. 
Eventually, a drops to zero at T = 0. This agrees with the result in Fig. 1 of [53] in 
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3+1 dimensions. It turns out that our model predicts an "insulator" at zero temperature. 
However, this is not a conventional insulator. In general, the insulator is gapped, i.e. 
the conductivity is effectively zero when the temperature scale is below the energy gap. 
However, the smooth behavior of the conductivity reducing to zero as T/fi — >■ in the top 
right plot of Fig. \5\ indicates that our model describes an exotic ungapped insulator at zero 
temperature. Since in our model, a depend only on the charged black brane background 
but not on 6, the zero temperature behavior has nothing to do with the break down 
of the pseudo scalar probe limit. [27] demonstrates that the scalar sector undergoes a 
smooth phase transition at T = 0, and according to our Fig. it turns out such phase 
transition related to the black brane background is also smooth and continuous. Fig. 
also displays different low temperature behavior of a compared to that in [28], in which 
cr is temperature independent in 2+1 dimensions. This is because our g vv has a second 
order zero at zero temperature, while the analysis in [2H] only applies to g vv with a first 
order zero. 

On the other hand, with p = 2Q, at high temperature a — > 1 and -A^ -+ This is 
different from the results obtained via linear response formalism in [29] where the electric 
conductivity in 2+1 dimensions is temperature independent, because the gauge fields in 
their analysis are in the probe limit, while our background is charged and already includes 
the backreaction. Our result match theirs at Q — > 0. 

4.2.2 Parity-odd part 

The equations ( l3Tj) and (132]) implies that ji(r) and ctj(r) take the form 

ji = jf ) + A 2 J f ) + 0(A 3 ), (43) 
ai = a f) + A 2 af ) + 0(A 3 ), 

The full vector mode equations (IBTi) and (j82j) in the Appendix show that, in the probe 
limit, the leading order 0(A°) equations satisfied by jf \ af' as in fl3T|) . (T32l describe 
only parity-even dynamics,. The subleading 0(A 2 ) equations satisfied by jf^ and af^ are 
composed of the parity-even and parity-odd parts. The former is regarded as the higher 
order correction to the leading order parity-even physics. The latter, however, gives rise 
to the leading order parity-odd effect in the vector mode. 

As a result, j\ 2 \ af^ can further be decomposed into jf^ = jf^ even ^ + j^dd) an< ^ = 
a f(iven) + a f{oddy wnere tne subscripts (ewen ) and {odd) label the parity. In order to study 
the parity violating hydrodynamics on the boundary, we have to solve the 0(\ 2 )( dd) bulk 
equations. We will denote j^ odd y a f( odd ) by ji, a-P'" 1 in the following. 

The two 0(A 2 )( o ^) equations of motion arising from -E^ 1 ' 2 ' 1 and r 2 M^ 1,2 \ after the 
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constraint (TH3T) is applied, are 

1 d 

2r 2 dr 



1 d (r%®') - lA%®> 



< ,;, djQ S Q + < i j djM 



1 1 d 



3M 4r 2 dr 
Si 

\r~ t a i := 

3MJ4r 2 dr K J ' 2r 2 



rj'6') - S 



M 



2 Q \ 1 ^ < r A f'0') 



d 



AQ 

-4 djQ - — g'^ij djM + 2 g'^j Ej 



S 2 



(44) 



(45) 



where -§\ and S 2 denote the sources of these differential equations, and 



Sm 
So 



2r 6 f 
1 



dM 4 

, de 



1 iff)' 



dM 4r 2 



r 4 # 



9 M {r) 



2r 6f> v ■> J qq 4 

9' Q 86 

r r 2 dQ ' 

& 3M d6 

r r 2 dM' 



1 (r 2 /')" 5 "' 1 



r 4 # 



dM 



dQ 4r 2 V y ' 



(46) 



3Mr 



- 0'. 



It is clear from above that the functions g<g, gu-, gE are no more dominant than 0(^33) 
as r 00. 

_____ /2"\ 

Integrating fT4"5]) from r# to r on both sides and then substituting d\ into (jSJ), one 



arrives at the differential equation for ji 



(2). 



d f 4 f2 d ( Ji i2 \r) \\ 
dr\ J dr\ f Ji 



4g 2 .(2) 



2Q 



jr(r H ) + fS 1 + ^ / s 



(47) 



in 



The solution is obtained by integrating the above equation from 00 to r on both sides 
twice, 



J^\r) = f 



s 4 f 2 



4Q 2 



(2) 

Then we can write down the expression of jf on the horizon: 



fs 1+ 2 4 r s 



111 



3M 



in 



fSi 



2Q r 



III 



(48) 



(49) 



This expression is consistent with the regularity condition. On the other hand, careful 
analysis shows that ji behaves asymptotically as 

h ~, 1 



^ (2) (^oc)^ + 0(^) + 
where • • • are terms of higher order than 0( ^+4 ), and 

h = -Q 2 Ji®(r H ) . 
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(50) 



(51) 



Since we limit ourself to < A_ < |, the leading order in floTJ]) is r 4 . Moreover, (15171) is 
consistent with the boundary conditions stated below (136]) . 

Then, we can evaluate the large r behavior of a\ by integrating (JHJ) on both side 
from oo to some large r as in (IBUj) . 

af ) (r->oo)~-~6 4 + .... (52) 
r Q 

The coefficients of 1/r in a\ (r — > oo) is interpreted as the expectation value correspond- 

— (2) 

ing to a\ on the boundary. Indeed, this is given by 



[Ji) = lim 



= lim yf^F ri ® = |- 6 4 

r— ><x> (J 

{odd) 



= K E eyEj + Kq eijdjQ + K M eijdjM . (53) 
After some algebra, one obtains 

Ke - ^{-mT^-smJ^ ^'mL (r/) \ 
R _ 4g K/'| 4g 2 r- to i r% 2nV \ 

Km " 3M \-^AT 9MlrH ~ 3M L L rfSM ~12ML ^ f ] * J ' 

in which we have carried out integration by part for the term J rH f(r 4 f'9')'. 

Let's compare the expression (I53|) with the parity violating part of the constitute 
equation of the boundary current J* in (jSJ): 

•/' = ?ey /'-■ + k (,,V, I + O'-'V//- . (55) 

Substituting in T, /i from the bulk theory in (|16p . one can identify the anomalous Hall 
conductivity a, anomalous thermal Hall conductivity k, and anomalous heat Hall con- 
ductivity £ in terms of Ke, Kq and Km- 



a = K E , 

1 3Mr H - 4Q 2 
~4nr 2 H 6Mr H - 4Q 2 



(2QK M + 3Mr hKq) , (56) 



« = 3Mr„-4Q^ 3M/?1 ' + 2 ^)- 

Since fj^, A'q and X M contain integrals of functions of 9 and 9', their values are 
subject to the choice of the boundary condition for 9. In the following we present the 
numerical result of a, k, and £ according to the two different types of boundary conditions 
we choose: 
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(i) Asymptotically sourced 9: 

As the source of 9 on the asymptotic boundary is turned on, the numerical results 
summarized in Fig. [6] ~ [8] shows that: a is positive and finite throughout who range of 
T/fi; is a l so positive and increases as T/fi reduces from large value before it starts to 
decrease below some small Tj\i. It eventually vanishes at T = as predicted by fl56|) ; £ 
is positive, increases monotonically as T j \x reduces, and diverges at T j \i = 0. 

(ii) Asymptotically sourceless 9: 

For the case where 9 is asymptotically sourceless and the boundary is parity sponta- 
neously broken by the nonzero vev of the dual pseudo scalar operator, Fig. I9KTTT1 display 
the numerical results of the three parity violating coefficients from the vector mode: a is 
nonvanishing and positive below the critical temperature T c , with a critical exponent 0.5 
near T c . -A^ and £ both diverge at T = T c and T = 0, both with critical exponents -1.5 
close to T = T c . These critical exponents, however, cannot be derived from fl56l) . as we 
don't have an analytical expression of 9. 



4.3 Tensor mode 



For the tensor mode, there is one dynamical equation due to HE X ]J . In the backreacted 
case it reads 



d_ 

dr 



lr 4 / d 



- r(d x (3 y + d y (3 Xj 



A d .r^fff 



Adr y H 2 



0. 



(57) 



4.3.1 Parity-even part 

After taking the probe limit 9 — >■ A# and p —> \<p, the 0(A°) part of the above equation 
becomes 

r(d x p y + d v p x ), (5? 



d 
dr 



2 5 dr xy 



in which the source only contains the shear mode. This equation can be analytically 
solved by integrating on both sides twice, with first integration being indefinite and the 
second one from oo to r: 

d x (3 y + d y (3 x f r 2a 3 



a 



(o) 
■ni 



r 2 f 



r 4f 



(59) 



The asymptotic behavior of a X y is 



a { xy \r^ oo) ~ 



d x Py + d y p x _^ 203 + 



(60) 



which is consistent with the sourceless boundary condition that the O(r ) coefficient 
vanishes. The 0(r~ 3 ) coefficient gives rise to the expectation value, 

3 2 



/ T (o) 



16/tGat 3 



«3, 



(61) 
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Figure 6: a versus T/n for various 9 mass m 2 , under the sourced boundary condition for 
9. The right figure is for small T / /i while the left one is for a wider range of T '/ ' \x. 



Mass square nr 




Figure 7: -jfh versus T/^u for various m 2 , under the sourced boundary condition for 9. 
The right figure is for small T/fi while the left one is for a wider range of T/fi. 




T 
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Figure 8: £ versus T/n for various m 2 , under the sourced boundary condition for 9. The 
right figure is for small T / /iwhile the left one is for a wider range of T/fi . 
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Figure 9: a versus T/fi with linear scale (on the left) and log scale (on the right) as 9 
is asymptotically sourceless. The left figure shows a vanishes at and above the critical 
temperature which depends on m 2 , and the right one indicates the critical exponent is 
0.5. 




M — 




Figure 10: ^ versus T//i with linear scale (on the left) and log scale (on the right) 
as 6 is asymptotically sourceless. The left figure shows that ^ diverges at the critical 
temperature, with a critical exponent —1.5 as being indicated by the right one. 
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Figure 11: £ versus T/fi with linear scale (on the left) and log scale (on the right) as 9 is 
asymptotically sourceless. The left figure shows that £ diverges at the critical temperature, 
with a critical exponent —1.5 as being indicated by the right one. 
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from which the shear viscosity r\ can be extracted. The constant is determined by the 
regularity condition on the horizon for a, 



(o) 



dr ^ v r 3 H f'(r H ) 



and CT: 



.2 



'r H f(r H ) rjj_ 
2 

Then one immediately obtains the shear viscosity 



* = Tfi^- (64) 

1D7TG7V 

r 2 

With the entropy density s = ^q- : the rj/s ratio match the universal bound, 

- = < 65 > 

4.3.2 Parity-odd part 

Next, we proceed to 0(A 2 ) order of the tensor mode equation (157|) . For this purpose we 
consider the a xy ansatz up to 0(A 2 ), 



a 



xy ^a® + \ 2 a®, (66) 
and substitute to (157|) . Then it is easy to see that the 0(A 2 ) equation is composed of 

(2) 

parity even and parity odd part, which implies that a xy can further be decomposed into 
two independent modes 

o§ = o$(r)(d m Py + d y /3 x ) + a®{r){d x p x - d y /3 y ), (67) 

where the shear mode is parity even and regarded as the 0(X 2 ) correction of the shear 
viscosity, while the Hall mode (proportional to d x (3 x — d y (3 y ) is the leading order parity odd 
effect. In the following we focus on the parity violating mode, a xy \ dd) = oixy (?*) (d x (3 x — 
d y P y ) and extract the Hall viscosity from bulk gravity. 
The 0(A 2 ) parity odd equation is 



d_ 

dr 



2 J dr xy 



\±ym ( 6 8) 



Integrating both sides twice, with first integration being indefinite and the second one 
from oo to r, one can write down the analytic solution for a 



xy j 



?& r 2d 



< w =>w-Lw (69) 
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Figure 12: tja/s versus T/fi for asymptotically sourced 9. The right figure is for small 
T/fi, while the left one for a wider range of T / fi. 



The asymptotic behavior is 



2«3 

3r 3 



+ 0( 



(70) 



This is consistent with the normalizable boundary condition. According to (1611) . the Hall 
viscosity r/A is encoded in a 3 which is determined by demanding (IBTIj) satisfy the regularity 
condition on the horizon, 



— a &) (r H ) 



1 1 d 
1 lr T f'~dr 



(71) 



Then a 3 is fixed to 



«3 



± 
4/ 



-(r'f'9'Y 



+ 



'7/ 



3M 



e'(r H ). 



(72) 



Such that 
Va 

VA 

s 



r 4 H f(r H )e'(r H ) 

1 dV 
' 8^~dJ rn = 



h 



'7/ 



327tGat 



(73) 



'7/ 



rn 



In the above expressions we have applied the regularity condition on the horizon for 9 in 
(TT8l) . It is clear from fT73l that the Hall viscosity is controlled by the value of 9 on the 
horizon as the mass square of the pseudo scalar. Since 9{th) is affected by the choice of 
the asymptotic boundary condition, so is Hall viscosity: 
(i) Asymptotically sourced 9: 

Fig. [121 displays tja/s as a function of T/fi, as 9 is subject to the condition that the 
asymptotic source on the boundary is turned on. The numerical analysis shows that, 
as the temperature reduces, rj A /s increases monotonically until the temperature is below 
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Figure 13: tja/s versus T/fi with linear scale (on the left) and log scale (on the right) 
for various 9 mass m 2 , as 9 is asymptotically sourceless. The left figure shows the Hall 
viscosity vanishes at and above some critical temperature, and the right figure indicates 
the critical exponent equal to 0.5. 



some small value, and then tja/s decreases along with the temperature and eventually 
reaches some finite value at T = 0. 
(ii) Asymptotically sourceless 9: 

Fig. [TBI shows the numerical results of 7]a/s versus T/fi when 9 is sourceless on the 
boundary. In this case, the Hall viscosity acquires non-vanishing value only below the 
critical temperature T c depending on m 2 . The critical exponent near T c is 0.5, as expected 
from the analysis of 9. 

We close this subsection by recalling the analytic solution of a xy for the general case 
with backreaction studied in ITBI: 



2H(t)dt 



a xy(r) = / ,A f u\ / dz 



t 4 f(t) 



I'll 



z(d x (3 y + d y (5 x ) + ^±(J^f)(d a /3 x - d yf 3 : 



(74) 



The above results of the shear and Hall viscosities are consistent with those derived from 
(1741) in [lj. 



5 Discussion and conclusion 

We have constructed a 3+1 dimensional holographic model with a charged black brane 
and neutral pseudo scalar field 9, dual to the parity violating hydrodynamic system in 2+1 
dimensions. We have extracted all the transport coefficients from the scalar, vector, and 
tensor modes of bulk perturbations, and reconstructed the hydrodynamics to first order in 
derivatives in the probe limit of 9. In our model, the parity violating effect is contributed 
from the neutral pseudo scalar 9, coupled to the gravitational and electrodynamical Chern- 
Simons terms. 
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The potential term in (jTJ allows the parity of the system to be broken spontaneously 
by the non-trivial vev of the boundary operator dual to 9 developed below T c , such that all 
parity transport coefficients are the order parameter for the second order phase transition 
and are switched on below T c . They exhibit certain critical exponent behavior very close to 
T c . Alternatively, one can choose the sourced boundary condition at asymptotic infinity, 
such that the parity and the conformal symmetry of the system is broken by this external 
source. The finite temperature behavior of all the transport coefficients arising in these 
two mechanisms are studied and presented. This paper generalize our previous work in 

One comment we would like to remind the readers is, in the case where the parity 
is spontaneously broken by 6, our result doesn't violate the Coleman-Mermin- Wagner 
theorem which states that in 2+1 dimensions the global continuous symmetry cannot 
be spontaneous broken, otherwise the correlation of the massless Goldstone boson would 
have infrared divergence. Since the spontaneous broken parity in our model is a discrete 
symmetry, no Goldstone boson is created. Moreover, our result shows that a and t}a/ s 
go to zero at T c with the same mean field exponent as the order parameter, (T c — T) a5 . 
However, -A-^ and £ exhibit (T c — T)~ L5 behavior near T c , which is novel and quite different 
from the usual systems. It would be interesting if these properties can be realized in the 
simple field theory models. 

Our bulk model has a U(l) gauge symmetry. It is interesting to generalize our analysis 
to the non-Abelian case, and construct dual hydrodynamics with non-Abelian gauge sym- 
metry. In 2+1 dimensions, the Hall conductivity can arise from a bare Chern-Simons term 
in the underlying theory, or from the induced one due to the quantum loop correction. It 
is known that the coupling constant of the non-Abelian CS term at T = has to be quan- 
tized in order to maintain the gauge symmetry. But in the finite temperature case, the 
coupling of the non-Abelian induced CS term is temperature dependent, and the gauge 
symmetry cannot be preserved for arbitrary temperature, unless all orders of quantum 
corrections are taken into account [36J. Therefore, non-perturbative methods are impor- 
tant in computing the Hall conductivity. Since we have used the gravity-hydrodynamics 
correspondence as an non-perturbative approach to derive the Hall conductivity, our anal- 
ysis may shed new lights on the non-perturbative calculation of the Hall conductivity for 
non-Abelian gauge theory. 

Very recently, [37] constructs a holographic model for 2+1 dimensional fluids with 
background vorticity Q and magnetic field B. They obtain the classical Hall conductivity 
inversely proportional to B, as expected. They also obtained the Hall viscosity inversely 
proportional to Q. This is different from our result, because in our model, we don't 
have non-vanishing external magnetic field and vorticity at the background level, i.e. in 
equilibrium; instead they appear as the first order perturbations. Therefore our Hall 
conductivity is "anomalous", and our Hall viscosity does not depend on vorticity. Their 
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result are derived by the Kubo's formula, and it is interesting to realize their result via 
gravity- hydrodynamics correspondence. 

Finally, We would like to point out that, in 2+1 dimensional (or more generally, in 
odd dimensional) quantum field theory, there exists the parity anomaly [3H] analogous 
to the even dimensional axial anomaly. As the case that the gauge anomaly induces 
new transport coefficients in the hydrodynamics [8], we expect to see the effect of this 
odd dimensional parity anomaly in the hydrodynamic transport too. We are studying 
the holographic realization of this effect in detail, and will present our result in the near 
future. 
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A Calculation of the first order scalar mode 

The equations of motion of the our bulk theory are given in ([8]). The general ansatz of 
black brane background with backreaction from the pseudo scalar field is ([9]). The stan- 
dard gravity-hydrodynamics approach requires to perturb the bulk system by promoting 
u^, M, Q, A e * 1 to the slow-varying functions of the boundary coordinates x^. Under the 
derivative expansions of the ansatz as in ( |T3|) without taking any probe limit, the 0(e°) 
ansatz satisfies the 0(e°) equations ([8]). The 0(e) part of the ansatz satisfies the 0(e) 
equations of motion which are given in this and the next section. 

For the scalar modes, there four dynamical equations at 0(e), arising from Erl\ 
H(r)M v ^\ r -^Er l r + E^) and KG^\ By applying the zeroth order equations of motion 
which also imply their derivatives by Q and M vanish, these four first order equations 
take the form 





(75) 



(76) 




<P + 



(77) 
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where ' denotes the derivative by r. There are also two constraint equations due to 
g rr E^r + g vr Evv and M r ^ which describe the relations between d v M, 8 V Q and difif 



r 8f 2rf B H + *MJ 



H 2 8M H 3 8M 2H 2 dM 



2H 2 



difr , (79) 



1 dH 1 dA'\ I \ A 1 

W8M A '-Wm) d " M + M -* Q ) d " Q = W m ■ (80) 



To proceed, we go into the probe limit 

9 -> X9, <f -»■ \<p, A -)• 
Then it follows the solutions in Section 14.11 



B Calculation of the first order vector mode 



For the 0(e) vector mode perturbations, there are two dynamical equations arising from 

E ri(l) and r 2 M i(l). 



1 d 

2r 2 dr 



4 • / 

H 



l_d_ ( r 2 fa< 
Hdr[ H 



+ a 

2H 1 



+ 



-r 2 d ( 1 8H\ 1 , d9 
~2~dr (rmdMj + 2 dM 



x (_H_d_ fr 7 f 2 



<9# 



1 d 



diM- 
r 2 f 



4Hdr\ H dM 



M - 
89' 



Q diQ 



1 d fr 4 9' 8f 



8M 

1 d fr 4 f'9'8H 
^~dr \ H 3 8M ) ^ 4r 2 dr \ H 2 
1 d fr 4 f'9' 
4r 2 c/r 



dj 8jM 



+A |m ->• q\ eij 8jQ + A 
r 2 A' , 



if 2 



# 2 



89 



H \8M 



A' - 



^9'\ } eij 8 3 m + \[m^ q) ey a,.g 



Off 

\— e ..p ext 
H m 



(82) 



28 



For both equations, the first line on the right hand side is parity-even, and the rest is 
parity-odd. After taking the probe limit, the leading order (0(A°)) terms contain only 
parity-even part, while the subleading (0(A 2 )) terms contain both parity-even and parity- 
odd contribution. 

There is one constraint equation from + E$: 

{ - — — (r^\ + —— (r 2 f) °* + - W 

1 2H dr V dM ) 2H 2 dr K J ' dM 2H dM 2H dM J 1 

+ |m Q JftQ = ^(r 2 f - AA') dvh - ^F:? . (83) 

By taking the probe limit and making use of /, H, A solutions, one can continue to 
the analysis in Section 14.21 The leading order gives rise to the parity-even transport 
coefficients a. The parity-even part of 0(X 2 ) equations are higher order corrections to 
the leading order hydrodynamics. The parity-even part of 0(A 2 ) gives rise to the parity 
violating coefficients. 
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